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We study one-particle spectra and a two-particle correlation function in the 130 GeV/nucleon 
Au+Au collisions at RHIC by making use of a hydrodynamical model. We calculate the one-particle 
hadronic spectra and present the first analysis of Bose-Einstein correlation functions based on the 
numerical solution of the hydrodynamical equations which takes both longitudinal and transverse 
expansion into account appropriately. The hydrodynamical model provides excellent agreement 
with the experimental data in the pseudorapidity and the transverse momentum spectra of charged 
hadrons, the rapidity dependence of anti-proton to proton ratio, and almost consistent result for 
the pion Bose-Einstein correlation functions. Our numerical solution with simple freeze-out picture 
, suggests the formation of the quark-gluon plasma with large volume and low net-baryon density. 



PACS numbers: 24.10.Nz, 12.38.Mh, 25.75.Gz 



Relativistic heavy ion collisions are very attracting problems which provide us the nature of hot and dense hadronic 
matter jjj. Creation of a new state of the matter, the quark-gluon plasma (QGP), and many kinds of new phenomena 
\ are expected to be found in the Relativistic Heavy Ion Collider (RHIC) experiments at BNL of which the collision 
energy is much higher than any other accelerator. However, the complicated processes during the many-body interac- 
^\ [ tions and multiparticle productions are quite hard to catch clear. Therefore, a simple phenomenological description is 
. indispensable for the better understanding of the phenomena. The aims of this paper are, based on a hydrodynamical 
model, to draw a simple and clear picture of the space-time evolution of the hot and dense matter produced in the 
high energy heavy ion collisions at RHIC and to give a possible explanation for the recent experimental results. 

We use a (3+l)-dimensional hydrodynamical model || to describe the space-time evolution assuming the local 
thermal and chemical equilibrium. Several authors have already discussed RHIC results based on hydrodynamical 



models. Kolb et al. |3] discussed anisotropic flow by making use of a (2+l)-dimcnsional hydrodynamic model in 
which Bjorken's ansatz [ El was used for the beam direction. Hence, their discussion is limited only in the midrapidity 
region. Zschieschc ct al.| p] discussed HBT radii based on a hydrodynamical model with use of the Bjorken's scaling 
solution in the collision direction. One of the authors (T.H.) has already reproduced the both the pseudorapidity and 
the transverse momentum spectra of hadrons by using a full (3+l)-dimensional hydrodynamical model in Ref. [p], 
where main theme of the analysis is also anisotropic flow. In this paper, we focus our discussion on central collisions 
, by assuming the cylindrical symmetry of the system. We calculate the one-particle hadronic spectra and present the 
k>( ' nrs ^ analysis of Bose-Einstein correlation functions based on the numerical solution of the hydrodynamical equations 
rS | which takes both longitudinal and transverse expansion into account appropriately 

The hydrodynamical equations are given as d tJ T^ l '{x) = with the baryon number conservation law d^n^x) = 0. 
We numerically solve these coupled equations for the perfect fluid by the method described in Ref. . Our numerical 
solution keeps entropy, energy and net baryon number conserved within 5% of accuracy throughout the calculation 
with the time step 5t — 0.01 fm/c. As for an equation of state (EOS), we adopt a bag model EOS in which phase 
transition of first order takes place 0. The QGP phase is a free gas with a bag constant B. The gas consists of 
massless quarks of three flavor and gluons. The hadronic phase is a free resonance gas with an excluded volume 
correction in which all resonances up to 2 GeV/c 2 of mass are included. These two phases are connected by the 
condition of pressure continuity 0. Putting the initial time as t — 1.0 fm/c, we parameterize the initial energy 
density distribution and net baryon number distribution as simple gaussian forms (Fig. |l|) Q. 1 The parameters in the 
model should be chosen so that calculated single-particle spectra reproduce the experimental results of Au+Au central 
collisions in the 130 GeV/nucleon at RHIC. The parameter set is summarized in Table. |. Note that we use Bjorken's 
solution 0] 1l = 77 only as an initial condition for the longitudinal flow velocity. Initial transverse flow is simply 
neglected. Once the freeze-out hypersurface is fixed, one can calculate the single particle spectra via Cooper-Frye 
formula m. We assume that the freeze-out process occurs a specific temperature Tf. This condition corresponds to 
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a constant energy density in the present calculation because of small baryonic chemical potential. By virtue of the 
Lagrangian hydrodynamics, we expect that contributions of the time- like hypersurface are small and the space- like 
hypersurface dominates the particle emission at freeze-out; we simply put k^da^ ~ k T da T . 2 We take hadrons from 
decay of resonances into account as well as directly emitted particles from the freeze-out hypersurface. We include 
decay processes p — > 2n, lo — > 37r, rj — > 37r, if* — * 7rii", and A — > 7V7T [ffol [n]. The resonances are also assumed to be 
emitted from the freeze-out hypersurface. Figure ^ shows pseudorapidity distribution of charged hadrons. Figure || 
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FIG. 2: Pseudorapidity %, distribution of charged hadrons. Solid line shows our result (n, K,p) including resonance contribution. 
Dotted line denotes contribution of the directly emitted particles from the freeze-out hypersurface. Closed circles are preliminary 
result from the PHOBOS Collaboration Q. 
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2 It is well known that the Cooper-Frye formula has an ambiguity in the treatment of the time-like hypersurface 
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FIG. 3: Transverse momentum spectrum of negatively charged hadrons. As in Fig. |2], solid line and dotted line show total 
number of particles and directly emitted particles from the freeze-out hypersurface, respectively. Closed circles are data from 
the STAR Collaboration 01. 
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FIG. 4: Transverse mass spectra of negatively charged hadrons. Solid line, dotted line and dashed line denote ir~ , K~ and p 
yield of our result. K~ and p spectra are scaled down by factor 0.1 and 0.01, respectively. Closed circles, open triangles and 
closed triangles are preliminary data from the PHENIX Collaboration |23|. 
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FIG. 5: Rapidity dependence of anti-proton to proton ratio. Experimental data are taken from the BRAHMS Collaboration 
& 

shows transverse momentum spectrum of negatively charged hadrons. In both figures, our results well reproduce the 
experimental data. We display the transverse mass spectra of identified negatively charged hadrons {~k~ , K~ , p) in 
Fig. [I| All slopes of the spectra are well reproduced by our calculation. For K~ and p, our results fail to reproduce 
the particle numbers but slopes agree well. In Fig. H, our result shows excellent agreement with the experimental 
data of anti-proton to proton ratio as a function of rapidity. However, the absolute numbers of anti-protons was 
not enough (Fig. |^). Though we here assume the same freeze-out condition for all particle species, the discrepancy 
may indicate the more complicated mechanism. From the parameter set in Table. [|, we can see that the energy 
density is sufficient for the QGP production. However, the energy density, e m ax = 6.0 GeV/fm 3 which corresponds 
to the temperature of 229 MeV, is only 5% higher than SPS fl^| . Large collision energy at RHIC leads to very large 
volume of the hot matter. Longitudinal extension of the hot matter, a ri + rjo = 2.47, is about 2.3 times larger than 
the one of SPS 0|. Initial energy density itself depends strongly on the initial time to which corresponds to the 
spatial size in the longitudinal direction. Hence, much higher energy density should be obtained if we assume an 
earlier thermalization time. Furthermore, we do not include the thickness of incident nuclei but use almost flat profile 
with gaussian smearing for the transverse direction (see Fig. |l|). Therefore, maximum energy density of our model 
is smaller than other models which take the thickness into account || EL n. The average energy density at rj = 
transverse plane in our model is 3.9 GeV/fm 3 . Total energy of the fluid of our model corresponds to 25290 GeV, 
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TABLE I: Parameter set for the Au+Au collisions. 



Maximum initial energy density Eq 6.0 GeV/fm 3 

Maximum initial net baryon density tibo 0.125 fm -3 

Longitudinal gaussian width a v of initial 1.47 
energy density 

Longitudinal extension r/o of the flat region 1.0 
in the initial energy density 

Longitudinal gaussian width <td of the ini- 1.4 
tial net baryon density 

Space-time rapidity 7/d at maximum of the 3.0 
initial net baryon distribution 

Gaussian smearing parameter a T of the 1.0 fm 

transverse profile 

Freeze-out temperature Tf 125 MeV 



TABLE II: Output. 



Net baryon number 131 

Mean chemical potential at freeze-out (/jb) 76.1 MeV 

Mean transverse flow velocity (vt) of the 0.509c 
fluid at \r)\ < 0.1 

Lifetime of the QGP phase tqgp 2.92 fm/c 

Lifetime of the mixed phase tmix 12.61 fm/c 

Total lifetime of the fluid thad 18.94 fm/c 



99 % of total collision energy. We display the outputs from the fluid in Table. ||. In the present model, net baryon 
number is much smaller than total baryon number of incident nuclei, (/^b) means average of the chemical potential 
on the whole freeze-out hypersurface. (v-p) is the average transverse velocity of the fluid element in |?7 < 0.1 at the 
freeze-out. "Lifetimes" of the each phase are also shown in Table. |J. Space-time evolution of the fluid is displayed in 
Fig. 0. The large area between T = 160 MeV and T = 158 MeV means large space-time volume of the mixed phase. 
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FIG. 6: Temperature contour plot on r-T plane. Left and right figures stand for r\ = and r\ — 3, respectively. 



The two-particle correlation function for chaotic source is calculated through 



\I( g n K^)\ 2 



5 



1.4 
1.3 
1.2 
1.1 
1 

0.9 



Hydro 
STAR 



0.05 

9ob [GeV/r 
1.4 
1.3 
1 1.2 

1 

0.9 



0.1 



1.4 
1.3 
1.2 
1.1 
1 

0.9 



Hydro 
STAR 





Hydro - 




STAR 














0.05 0.1 

?long [GeV/c] 



0.05 0.1 

q siic [GeV/c] 



0.15 



0.15 



FIG. 7: Two-particle Bose-Einstein correlation functions for n~ . Upper-left, upper-right and lower figures show outward, 
sideward and longitudinal correlation functions, respectively. In each figure, correlation function is integrated with respect to 
other two components from to 35 MeV/c. Experimental data (closed circles) are taken from the STAR Collaboration |Ba|. 
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FIG. 8: Mr dependence of n~ HBT radii. Upper-left, upper-right, lower-left and lower-right figures correspond to outward, 
sideward, longitudinal HBT radii and ratio of outward HBT radii to sideward HBT radii, respectively. Experimental data 
(closed circles) are taken from the STAR Collaboration [BSJ. 



where K» = (fcf + k%)/2 t cf = fef - k%, respectively [§]§ 0. We put 

I(<f,K») = J K ti d^(xWf(k 1 ,x)f(k 2 ,x)e i ^ xU , (2) 

so that 1(0, /c M ) reduces to the Cooper-Frye formula with f(k, x) being the Bose-Einstein distribution function. As 
a first trial, we neglect the contributions from resonance decay for simplicity. Figure shows the projected correla- 
tion functions for 7r~ with \Y\ < 0.5 and 0.125 < Kt < 0.225 GeV/c. In each correlation function, our results are 
corrected by a common A factor whose value is 0.6. The other origin of the reduction is integration with respect 
to other components of the relative momenta over the range < < 35 MeV/c. Despite the neglection of reso- 
nance contribution, 3 both outward and longitudinal correlation functions show good agreements with the experiment. 
Finally, we compare pair transverse mass Mx = y K\ + m 2 dependence of the HBT radii extracted through the 
Gaussian fit to the three-dimensional correlation function (Fig. ||). We obtain almost consistent sideward HBT radii 



3 It has been discussed that the hadronic cloud significantly affects the HBT radii |15 



6 



but obtained outward and longitudinal HBT radii are larger than experimental data. This tendency is common to 
the all hydrodynamical calculation However, the difference between our results and data, 1 fm in HBT radii, 

corresponds to only 5 MeV in the correlation function (Fig. |?]). As for i? out , a naive interpretation of the experimental 
result is that high Mt pions come from high temperature source and the experimental results indicate short lifetime 
of the high temperature region. On the other hand, if the freeze-out picture works well, i.e., most of particles are 
emitted from the 3-dimensional hypcrsurface, the experimental data suggest very short freeze-out duration or strong 
opaqueness of the source 0. Though our model naturally shows opaque property due to hydrodynamical flow |L8| , 
the freeze-out time duration, At = \J (t 2 ) — (t) 2 where (A(x)} = J d A xA(x)S(x, K)/ J d 4 xS(x,K) with S(x, K) being 
the source function of our model is about 5 ~ 7 fm/c which is longer than the SPS case jL2|, fL8| , This long time 
duration of our model leads to the large -R ou t/-Rsidc @- 

In summary, we present a hydrodynamical-model calculation for the 130 GeV/nucleon Au+Au collisions data from 
the RHIC experiments. The present numerical solution indicates that the produced quark-gluon plasma in the RHIC 
has much lower baryon density, slightly higher energy density, and several times larger extension in the longitudinal 
direction than in the SPS case, if we compare at the same initial time, ro = 1 fm/c. We present the first analysis of 
the pion Bose-Einstein correlation data at RHIC based on the hydrodynamical model which takes into account both 
longitudinal and transverse expansion appropriately. The result is partly consistent with the experimental result but 
obtained i? ut and i?i ong are slightly larger than the experimental data. In this work, we concentrate our discussion on 
the RHIC data. More detailed discussion including the comparison with the SPS data will be given in the forthcoming 
paper @. 
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